Disequilibrating azobenzenes by visible-light sensitization under confinement.

Photoisomerization of azobenzenes from their stable E isomer to the metastable Z state is the basis of numerous applications of these molecules. However, this reaction typically requires ultraviolet light, which limits applicability. In this study, we introduce disequilibration by sensitization under confinement (DESC), a supramolecular approach to induce the E-to-Z isomerization by using light of a desired color, including red. DESC relies on a combination of a macrocyclic host and a photosensitizer, which act together to selectively bind and sensitize E-azobenzenes for isomerization. The Z isomer lacks strong affinity for and is expelled from the host, which can then convert additional E-azobenzenes to the Z state. In this way, the host-photosensitizer complex converts photon energy into chemical energy in the form of out-of-equilibrium photostationary states, including ones that cannot be accessed through direct photoexcitation.

Altered Nucleotide Insertion Mechanisms of Disease-Associated TERT Variants.

Telomere biology disorders (TBDs) are a spectrum of diseases that arise from mutations in genes responsible for maintaining telomere integrity. Human telomerase reverse transcriptase (hTERT) adds nucleotides to chromosome ends and is frequently mutated in individuals with TBDs. Previous studies have provided insight into how relative changes in hTERT activity can lead to pathological outcomes. However, the underlying mechanisms describing how disease-associated variants alter the physicochemical steps of nucleotide insertion remain poorly understood. To address this, we applied single-turnover kinetics and computer simulations to the Tribolium castaneum TERT (tcTERT) model system and characterized the nucleotide insertion mechanisms of six disease-associated variants. Each variant had distinct consequences on tcTERT's nucleotide insertion mechanism, including changes in nucleotide binding affinity, rates of catalysis, or ribonucleotide selectivity. Our computer simulations provide insight into how each variant disrupts active site organization, such as suboptimal positioning of active site residues, destabilization of the DNA 3' terminus, or changes in nucleotide sugar pucker. Collectively, this work provides a holistic characterization of the nucleotide insertion mechanisms for multiple disease-associated TERT variants and identifies additional functions of key active site residues during nucleotide insertion.

Future COVID19 surges prediction based on SARS-CoV-2 mutations surveillance.

COVID19 has aptly revealed that airborne viruses such as SARS-CoV-2 with the ability to rapidly mutate combined with high rates of transmission and fatality can cause a deadly worldwide pandemic in a matter of weeks (Plato et al., 2021). Apart from vaccines and post-infection treatment options, strategies for preparedness will be vital in responding to the current and future pandemics. Therefore, there is wide interest in approaches that allow predictions of increase in infections ('surges') before they occur. We describe here real-time genomic surveillance particularly based on mutation analysis, of viral proteins as a methodology for a priori determination of surge in number of infection cases. The full results are available for SARS-CoV-2 at http://pandemics.okstate.edu/covid19/, and are updated daily as new virus sequences become available. This approach is generic and will also be applicable to other pathogens.

The involvement of triplet states in the isomerization of retinaloids.

Rhodopsins form a family of photoreceptor proteins which utilize the retinal chromophore for light energy conversion. Upon light absorption the retinal chromophore undergoes a photoisomerization. This reaction involves a non-radiative relaxation through a conical intersection between the singlet excited state and the ground state. In this work we studied the possible involvement of triplet states in the photoisomerization of retinaloids using the extended multistate (XMS) version of CASPT2. To this end, truncated models of three retinaloids were considered: protonated Schiff base, deprotonated Schiff base and the aldehyde form. The optimized geometries of the reactant, the product and the conical intersection were connected by a linear interpolation of internal coordinates to describe the isomerization. The energetic position of the low-lying singlet and triplet states as well as their spin-orbit coupling matrix elements (SOCME) were calculated along the isomerization profile. The SOCME values peaked in vicinity of the conical intersection for all the retinaloids. Furthermore, the magnitude of SOCME is invariant to the number of double bonds in the model. The SOCME for the protonated Schiff base is negligible (1.5 cm-1) which renders the involvement of the triplet state as improbable. However, the largest SOCME value of 30 cm-1 was found for the aldehyde form, followed by 15 cm-1 for the deprotonated Schiff base.

Rhodopsin-bestrophin fusion proteins from unicellular algae form gigantic pentameric ion channels.

Many organisms sense light using rhodopsins, photoreceptive proteins containing a retinal chromophore. Here we report the discovery, structure and biophysical characterization of bestrhodopsins, a microbial rhodopsin subfamily from marine unicellular algae, in which one rhodopsin domain of eight transmembrane helices or, more often, two such domains in tandem, are C-terminally fused to a bestrophin channel. Cryo-EM analysis of a rhodopsin-rhodopsin-bestrophin fusion revealed that it forms a pentameric megacomplex (~700 kDa) with five rhodopsin pseudodimers surrounding the channel in the center. Bestrhodopsins are metastable and undergo photoconversion between red- and green-absorbing or green- and UVA-absorbing forms in the different variants. The retinal chromophore, in a unique binding pocket, photoisomerizes from all-trans to 11-cis form. Heterologously expressed bestrhodopsin behaves as a light-modulated anion channel.

Signal transduction in light-oxygen-voltage receptors lacking the active-site glutamine.

In nature as in biotechnology, light-oxygen-voltage photoreceptors perceive blue light to elicit spatiotemporally defined cellular responses. Photon absorption drives thioadduct formation between a conserved cysteine and the flavin chromophore. An equally conserved, proximal glutamine processes the resultant flavin protonation into downstream hydrogen-bond rearrangements. Here, we report that this glutamine, long deemed essential, is generally dispensable. In its absence, several light-oxygen-voltage receptors invariably retained productive, if often attenuated, signaling responses. Structures of a light-oxygen-voltage paradigm at around 1 Å resolution revealed highly similar light-induced conformational changes, irrespective of whether the glutamine is present. Naturally occurring, glutamine-deficient light-oxygen-voltage receptors likely serve as bona fide photoreceptors, as we showcase for a diguanylate cyclase. We propose that without the glutamine, water molecules transiently approach the chromophore and thus propagate flavin protonation downstream. Signaling without glutamine appears intrinsic to light-oxygen-voltage receptors, which pertains to biotechnological applications and suggests evolutionary descendance from redox-active flavoproteins.

Deciphering the Spectral Tuning Mechanism in Proteorhodopsin: The Dominant Role of Electrostatics Instead of Chromophore Geometry.

Proteorhodopsin (PR) is a photoactive proton pump found in marine bacteria. There are two phenotypes of PR exhibiting an environmental adaptation to the ocean's depth which tunes their maximum absorption: blue-absorbing proteorhodopsin (BPR) and green-absorbing proteorhodopsin (GPR). This blue/green color-shift is controlled by a glutamine to leucine substitution at position 105 which accounts for a 20 nm shift. Typically, spectral tuning in rhodopsins is rationalized by the external point charge model but the Q105L mutation is charge neutral. To study this tuning mechanism, we employed the hybrid QM/MM method with sampling from molecular dynamics. Our results reveal that the positive partial charge of glutamine near the C14 -C15 bond of retinal shortens the effective conjugation length of the chromophore compared to the leucine residue. The derived mechanism can be applied to explain the color regulation in other retinal proteins and can serve as a guideline for rational design of spectral shifts.

The impact of twisting on the intersystem crossing in acenes: an experimental and computational study.

Due to their unique excited state dynamics, acenes play a dominant role in optoelectronic and light-harvesting applications. Their optical and electronic properties are typically tailored by side-group engineering, which often result in distortion of the acene core from planarity. However, the effect of such distortion on their excited state dynamics is not clear. In this work, we investigate the effect of twisting on the photophysics of acenes, which are helically locked to a defined twist angle by tethers of different lengths. Ultrafast transient absorption and time resolved fluorescence show a clear dependence of the rate of intersystem crossing with twisting. This trend is explained using quantum chemical calculations, showing an increase of spin-orbit coupling (SOC). At much earlier times, structural reorganization in S1, including coherent vibrational wave packet motions, is reflected in transient spectral changes. As predicted by theory, decreasing the length of diagonal tether induces enhanced activity and frequency blue-shifting of a normal vibration consisting of anthracene twisting against restraint of the tethering chain. Overall, these results serve as design principles for tuning photophysical properties of acenes via controlled twisting of their aromatic core.

Bending versus Twisting Acenes - A Computational Study.

Polycyclic aromatic hydrocarbons (PAHs) are widely used in organic electronic devices. The electronic, magnetic, and optical properties of PAHs can be tuned by structural modifications to the aromatic backbone to introduce an inherent distortion from planarity, such as bending or twisting. However, it remains difficult to isolate and control the effects of such distortions. Here, we sought to understand how backbone twisting and bending affect the electronic properties of acenes, as models for larger PAHs. We found that, even when highly distorted from planarity (30° per ring), acenes maintain their aromatic character and π orbital delocalization with minor mixing of the σ and π orbitals. In addition, the energy gap between the HOMO and LUMO decreases with increasing twist, while the gap is hardly affected by bending, since the energy of both orbitals increase to a similar extent. For bent acenes in the triplet state, the spin becomes more localized with increasing bend, whereas twisting produces an evenly distributed spin delocalization. These findings can guide the synthesis of PAHs with tailored properties.

Frontiers in Multiscale Modeling of Photoreceptor Proteins.

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

NeoR, a near-infrared absorbing rhodopsin.

The Rhizoclosmatium globosum genome encodes three rhodopsin-guanylyl cyclases (RGCs), which are predicted to facilitate visual orientation of the fungal zoospores. Here, we show that RGC1 and RGC2 function as light-activated cyclases only upon heterodimerization with RGC3 (NeoR). RGC1/2 utilize conventional green or blue-light-sensitive rhodopsins (λmax = 550 and 480 nm, respectively), with short-lived signaling states, responsible for light-activation of the enzyme. The bistable NeoR is photoswitchable between a near-infrared-sensitive (NIR, λmax = 690 nm) highly fluorescent state (QF = 0.2) and a UV-sensitive non-fluorescent state, thereby modulating the activity by NIR pre-illumination. No other rhodopsin has been reported so far to be functional as a heterooligomer, or as having such a long wavelength absorption or high fluorescence yield. Site-specific mutagenesis and hybrid quantum mechanics/molecular mechanics simulations support the idea that the unusual photochemical properties result from the rigidity of the retinal chromophore and a unique counterion triad composed of two glutamic and one aspartic acids. These findings substantially expand our understanding of the natural potential and limitations of spectral tuning in rhodopsin photoreceptors.

Modern quantum chemistry with [Open]Molcas.

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Publisher Correction: Roaming-mediated ultrafast isomerization of geminal tri-bromides in the gas and liquid phases.

In the version of this Article originally published, in Fig. 5, the chemical formula Br•CC6H11 should have read Br•CH3C6H11.

Spectroscopic Properties of Lumiflavin: A Quantum Chemical Study.

In this work, the electronic structure and spectroscopic properties of lumiflavin are calculated using various quantum chemical methods. The excitation energies for ten singlet and triplet states as well as the analysis of the electron density difference are assessed using various wave function-based methods and density functionals. The relative order of singlet and triplet excited states is established on the basis of the coupled cluster method CC2. We find that at least seven singlet excited states are required to assign all peaks in the UV/Vis spectrum. In addition, we have studied the solvatochromic effect on the excitation energies and found differential effects except for the first bright excited state. Vibrational frequencies as well as IR, Raman and resonance Raman intensities are simulated and compared to their experimental counterparts. We have assigned peaks, assessed the effect of anharmonicity, and confirmed the previous assignments in case of the most intense transitions. Finally, we have studied the NMR shieldings and established the effect of the solvent polarity. The present study provides data for lumiflavin in the gas phase and in implicit solvent model that can be used as a reference for the protein-embedded flavin simulations and assignment of experimental spectra.

Retinal isomerization in bacteriorhodopsin captured by a femtosecond x-ray laser.

Ultrafast isomerization of retinal is the primary step in photoresponsive biological functions including vision in humans and ion transport across bacterial membranes. We used an x-ray laser to study the subpicosecond structural dynamics of retinal isomerization in the light-driven proton pump bacteriorhodopsin. A series of structural snapshots with near-atomic spatial resolution and temporal resolution in the femtosecond regime show how the excited all-trans retinal samples conformational states within the protein binding pocket before passing through a twisted geometry and emerging in the 13-cis conformation. Our findings suggest ultrafast collective motions of aspartic acid residues and functional water molecules in the proximity of the retinal Schiff base as a key facet of this stereoselective and efficient photochemical reaction.

A QM/MM study of the initial excited state dynamics of green-absorbing proteorhodopsin.

The primary photochemical reaction of the green-absorbing proteorhodopsin is studied by means of a hybrid quantum mechanics/molecular mechanics (QM/MM) approach. The simulations are based on a homology model derived from the blue-absorbing proteorhodopsin crystal structure. The geometry of retinal and the surrounding sidechains in the protein binding pocket were optimized using the QM/MM method. Starting from this geometry the isomerization was studied with a relaxed scan along the C13[double bond, length as m-dash]C14 dihedral. It revealed an "aborted bicycle pedal" mechanism of isomerization that was originally proposed by Warshel for bovine rhodopsin and bacteriorhodopsin. However, the isomerization involved the concerted rotation about C13[double bond, length as m-dash]C14 and C15[double bond, length as m-dash]N, with the latter being highly twisted but not isomerized. Further, the simulation showed an increased steric interaction between the hydrogen at the C14 of the isomerizing bond and the hydroxyl group at the neighbouring tyrosine 200. In addition, we have simulated a nonadiabatic trajectory which showed the timing of the isomerization. In the first 20 fs upon excitation the order of the conjugated double and single bonds is inverted, consecutively the C13[double bond, length as m-dash]C14 rotation is activated for 200 fs until the S1-S0 transition is detected. However, the isomerization is reverted due to the specific interaction with the tyrosine as observed along the relaxed scan calculation. Our simulations indicate that the retinal - tyrosine 200 interaction plays an important role in the outcome of the photoisomerization.

Excitation Energies of Canonical Nucleobases Computed by Multiconfigurational Perturbation Theories.

In this computational work, we assessed the performance of ab initio multireference (MR) methods for the calculation of vertical excitation energies of five nucleobases: adenine, guanine, cytosine, thymine and uracil. In total, we have studied 38 singlet and 30 triplet excited states. Where possible we used the multireference configuration interaction (MRCI) method as a reference for various flavors of multireference perturbation theory to second order. In particular, we have benchmarked CASPT2, NEVPT2 and XMCQDPT2. For CASPT2, we have analyzed the single-state, multistate (MS) and extended MS variants. In addition, we have assessed the effect of the ionization potential electron affinity (IPEA) shift. For NEVPT2, we have used the partially and the strongly contracted variants. Further, we have tested the commonly used RI-CC2, RI-ADC2 and EOM-CCSD methods. Generally, we observe the following trends for singlet excited states: NEVPT2 is the closest MR method to MRCISD+Q, closely followed by CASPT2 with the default IPEA shift. The same trend is observed for triplet states, although NEVPT2 and CASPT2-IPEA are getting closer. Interestingly, the n, π* singlet excited states were described more accurately than π, π* excited states, while for triplet states the trend is inverted except for NEVPT2. This work is an important benchmark for future photochemical investigations.

Direct photoisomerization of CH2I2vs. CHBr3 in the gas phase: a joint 50 fs experimental and multireference resonance-theoretical study.

Femtosecond transient absorption measurements powered by 40 fs laser pulses reveal that ultrafast isomerization takes place upon S1 excitation of both CH2I2 and CHBr3 in the gas phase. The photochemical conversion process is direct and intramolecular, i.e., it proceeds without caging media that have long been implicated in the photo-induced isomerization of polyhalogenated alkanes in condensed phases. Using multistate complete active space second order perturbation theory (MS-CASPT2) calculations, we investigate the structure of the photochemical reaction paths connecting the photoexcited species to their corresponding isomeric forms. Unconstrained minimum energy paths computed starting from the S1 Franck-Condon points lead to S1/S0 conical intersections, which directly connect the parent CHBr3 and CH2I2 molecules to their isomeric forms. Changes in the chemical bonding picture along the S1/S0 isomerization reaction path are described using multireference average coupled pair functional (MRACPF) calculations in conjunction with natural resonance theory (NRT) analysis. These calculations reveal a complex interplay between covalent, radical, ylidic, and ion-pair dominant resonance structures throughout the nonadiabatic photochemical isomerization processes described in this work.

Roaming-mediated ultrafast isomerization of geminal tri-bromides in the gas and liquid phases.

'Roaming' is a new and unusual class of reaction mechanism that has recently been discovered in unimolecular dissociation reactions of isolated molecules in the gas phase. It is characterized by frustrated bond cleavage, after which the two incipient fragments 'roam' on a flat region of the potential energy surface before reacting with one another. Here, we provide evidence that supports roaming in the liquid phase. We are now able to explain previous solution-phase experiments by comparing them with new ultrafast transient absorption data showing the photoisomerization of gas-phase CHBr3. We see that, upon S0-S1 excitation, gas-phase CHBr3 isomerizes within 100 fs into the BrHCBr-Br species, which is identical to what has been observed in solution. Similar sub-100 fs isomerization is now also observed for BBr3 and PBr3 in solution upon S1 excitation. Quantum chemical simulations of XBr3 (X = B, P or CH) suggest that photochemical reactivity in all three cases studied is governed by S1/S0 conical intersections and can best be described as occurring through roaming-mediated pathways.